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Abstract 

Tvo methods for modeling dynamic stall have 
been developed at United Aircraft. The a. A, B 
Method generates lift and pitching moments as 
functions of angle of attack and its first two 
time derivatives. The coefficients are derived 
from experimental data for oscillating airfoils. 
The Time Delay Method generates the coefficients 
from steady state airfoil characteristics and an 
associated time delay in stall beyond the steady 
state stall angle. Correlation with three types 
of test data shows that the a. A, B Method is 
somewhat better for use in predicting helicopter 
rotor response in forward flight. Correlation 
with lift and moment hysteresis loops generated 
for oscillating airfoils was good for both models. 
Correlation with test data in which flexibly 
mounted two-dimensional airfoils were oscillated 
to simulate the IP pitch variation of a helicopter 
rotor blade showed that both methods overpredicted 
the response, and neither gave a clear advantage. 
The a. A, B Method gave better correlation of 
torsional response of full scale rotors and re- 
mains the method in general use. The Time Delay 
Method has the potential to be applied more easily 
and probably can be improved by consideration of 
spanwise propagation of stall effects. 

Stall-related phenomena limit the operation- 
al capabilities of the helicopter. Power, blade 
stress , and control system loads can all increase 
substantially due to blade stall. To predict 
such phenomena unsteady aerodynamics in stall must 
be modeled in blace aeroelastic analyses . A num- 
ber of unsteady aerodynamic models have been 
developed. These include methods described in 
References 1 and 2. Reference 3 is a recent 
general survey article of rotor dynamic stall. 

The a. A, B Method and the Time Delay Method are 
two methods developed by United Aircraft. The 
a. A, B Method was developed to use airfoil test 
data obtained for a sinusoidally oscillating 
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two-dimensional model airfoil. The Time Delay 
Method was developed to provide an empirical method 
that would agree with the lift and pitching moment 
hysteresis char act eristics measured in oscillating 
airfoil tests for a number of airfoils and test 
conditions . 

Evaluation of unsteady aerodynamic modeling 
techniques generally proceeds from correlation with 
data obtained in two-dimensional oscillating air- 
foil tests to correlation of full scale rotor blade 
torsional response. Two-dimensional rigid airfoil 
results are compared on the basis of aerodynamic 
pitch damping and lift and pitching moment hys- 
teresis loops , and full scale correlation is judged 
on the basis of agreement in blade torsional mo- 
ments or control rod loads. Evaluation of an un- 
steady model on the basis of full scale torsional 
response is made difficult by uncertainties in 
three-dimensional rotor inflow and blade bending 
and plunging motion. Correlation of the lift and 
pitching moment time histories of rigidly driven 
airfoils, on the other hand, is not the best method 
of comparison because it does not treat blade dy- 
namic response to stall. As an intermediate ap- 
proach, model test data were obtained for a flex- 
ibly mounted model airfoil which was dynamically 
scaled to simulate the dynamics of the first tor- 
sional mode of a rotor blade. This paper summa- 
rizes unsteady aerodynamic modeling techniques and 
includes comparisons based on two-dimensional aero- 
dynamic pitch damping, lift and pitching moment 
hysteresis loops, two-dimensional flexured airfoil 
response, and full scale rotor blade torsional 
moments . 

Description of the Unsteady Models 
ot. A, B Method 

In the a. A, B method the aerodynamic moment 
is assumed to he a function of angle of attack and 
its first two time derivatives. Reference k demon- 
strated that unsteady normal force and moment data 
generated during sinusoidal airfoil tests and tabu- 
lated as functions of a, A = ba and B = b^a 

, u 0 

(where b is the airfoil semi-chord and U 0 is the 
free stream velocity) could he used to predict the 
aerodynamic response of an airfoil executing 
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a nonsinusoidal motion. In a limited number of 
flexured airfoil tests described in Reference 4, 
good correlation was achieved between measured 
and predicted airfoil dynamic response. The a, 

A, B lift and pitching moment data tabulations 
of Reference b were used in the calculation of 
torsional response for the dynamically scaled 
model airfoil. As applied in this investigation, 
two changes were made in the calculation. First, 
to consider the pitch axis of the model airfoil as 
a vari able, provision was made to include pitching 
moment due to chordwise offset of the aerodynamic 
center from the pitch axis: 

( a,A,B)=c m ( ot, A,B ) + ( Xp^~ .25 ( ot,A,B) 

7 c/4 

X PA 


break and then a loss of lift corresponding to flow 
separation. The relative phasing of the moment and 
lift breaks with angle of attack produces either 
positive or negative damping of the motion. 

To test the Time Delay hypothesis, harmonic 
data from Reference 5 were examined. It was noted 
that the onset of stall can occur before, with, or 
after maximum amplitude of the oscillation. In 
accordance with the Time Delay hypothesis, the 
spread between the static moment stall angle and 
the dynamic lift break was evaluated in terms of 
elapsed time nondimens ionali zed by free stream 
velocity and chord length, x# = At SE p(U o /c) • 

Typical results show that separation generally 
occurs when x# exceeds about 6 . 


The second change involved scaling the un- 
steady data tables to account for differences in 
wind tunnel characteristics. The steady state 
lift and moment data for the present test program 
differed from the corresponding steady data 
obtained in Reference 4 because the tests were 
conducted in an open jet wind tunnel and because 
the airfoil effective aspect ratio was much 
higher. The method of scaling used for these 
analyses required a shift in the entire data tabu- 
lation by constant values of angle of attack, un- 
steady lift coefficient, and unsteady moment co- 
efficient according to the following relations: 

c l ( “’ A ’ B) open jet =c l (a+ 6 a l’ A ’ B) TAB + 6c l 


and 


c (a,A,B) =c (a+6a ,A,B) m . + <5c 

m V* 'open jet m m* * TAB » 


The constants 60 ^, 6 cj_ and 6 c m were es- 

tablished for each airfoil and were equal to the 
amount of shift necessary to make the open jet 
steady state stall points in lift and moment 
match the steady state stall points of the 
airfoil of Reference 4. 


Time Delay Unsteady Model 

Wind tunnel airfoil dynamic response was 
also calculated with the Sikorsky Time Delay un- 
steady aerodynamic method. This formulation was 
developed empirically by generalizing the re- 
sults of a set of oscillating airfoil test pro- 
grams. It is intended to predict the unsteady 
aerodynamic characteristics of arbitrary airfoils. 
Its aim is to provide the blade designer with un- 
steady lift and pitching moment characteristics 
of various airfoils without conducting extensive 
oscillating airfoil tests. This model, based on 
a hypothesis of the physical separation process, 
does not depend on an assumed harmonic variation 
of angle of attack. The basic assumption is that 
there exists a maximum quasi-static angle of 
attack at which the pressure distribution and the 
boundary layer are in equilibrium. During in- 
creases in angle of attack beyond this static 
stall angle, there are finite time delays before 
a redistribution of pressure causes first a moment 


Dynamic pitching moment stall has been 
handled similarly. Test data showed, in general, 
that the dynamic moment break occurred before the 
lift break. This has been noted in Reference 6 
and attributed to the shedding of a vortex at the 
airfoil leading edge at the beginning of the 
separation process. Rearward movement of the vor- 
tex over the surface of the airfoil tends to main- 
tain lift, but drastically alters the pitching 
moment . 

To apply the Time Delay Model to a given air- 
foil, only static aerodynamic data are required. 
First, the airfoil static lift and pitching moment 
data are used to define the approximate variation 
in center of pressure between the static moment 
stall angle a and an angle of attack above 
which the center of pressure is assumed fixed. 
Secondly, an approximation is made to the c 
versus a curve for fully separated flow. The se- 
quence of events occurring during one stall-unstall 
cycle is detailed in Figure 1. Briefly stated, 
lift and pitching moment are determined from po- 
tential flow theory until the nondimensional time 
(which begins counting when the angle of attack 

exceeds the static moment stall angle) reaches x g . 
At this point the pressure distribution begins to 
change, leading to rearward movement of the center 
of pressure and loss of potential flow pitching 
moment. Later, when x^ = x*, the lift breaks from 

the static line and decreases gradually with time 
to the fully separated value, c.. (a). For 

dEP 

t 2 >t* the center of pressure coincides with 
C.P. SEp (a). At the point where a = 0, the rates 
at which approaches c i SE p( a ) C.P. approaches 

C.P.g Ep (a) (if it does not already equal C.P. gE p(a)) 

are increased. When a falls back below the quasi- 
static stall angle ot-p the center of pressure 
returns to the quarter chord, potential flow 
pitching moment effects are reintroduced and a 
second time parameter x is recorded to govern the 
rate at which c returns to c 
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Figure 1. Time Delay Unsteady Aerodynamic Model. 
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Figure 2. Correlation with NACA 0012 Lift and 
Pitching Moment Hysteresis Loops. 


Although additional correlation studies must 
be made to identify the effects of airfoil type on 
the time delay constant a and although refinements 
to the present model may be implemented, this 
rather ainple model represents well the essential 
features of the dynamic stall process. Correlation 
typical of that claimed for other empirical methods 
(Inferences 2 and 7) has been found with data from 
inferences 4 , 5, and 6. Only the a, A, B Method 
has produced better correlation (inference 4), but 
it Buffers from the requirement for extensive 
testing and data processing. Figure 2 compares 
the NACA 0012 unsteady lift and pitching moment 
hysteresis loops measured in Reference 4 with Time 
Daisy results. This correlation was achieved by 
setting the lift break time constant t* equal to 
4.0 instead of 6.0. Three-dimensional offsets 
encountered in this test apparently reduced the 
time interval between static stall and dynamic 
lift stall. Also shown are the hysteresis loops 
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Figure 3. Correlation of Two-Dimensional Aero- 
Dynamic Pitch Damping. 

predicted using the a, A, B Method. The a, A, B 
Method correlation is with the data from which the 
a. A, B coefficients were derived. 

In addition to predicting the exact form of 
lift and moment hysteresis loops, an unsteady- 
model should represent faithfully aerodynamic 
pitch: damping. Accordingly, the Time Delay 
Model was used to calculate two-dimensional aero- 
dynamic damping for the reduced frequency /me an 
angles of attack test points of Reference 9 - 
Sample results plotted versus airfoil mean in- 
cidence angle of attack are shown in Figure 3. 
Generally excellent correlation of measured and 
predicted damping is noted. 

Other correlation of the Time Delay Method 
with two-dimensional oscillating airfoil test 
data has been good. During development of the 
theory, correlation was carried out with forced 
oscillating airfoil data for a range of airfoils, 
frequencies of forced oscillation, Mach numbers, 
and angles of attack. Typical examples of the 
correlation obtained are shown in Figure 4 where 
measured and calculated hysteresis loops are 
shown for the V13006-7 airfoil. These test data 
taken from Reference 1 show the correlation with 
the Boeing Theory of Reference 1 as well. Corre- 
lation included hysteresis loops for different 
airfoils and covered a Mach number range from 0.2 
to 0.6. In all cases, the general character and 
magnitude of the hysteresis loops were well match- 
ed. In particular, the method provides the sharp 
drop in pitching moment that is often found when 
stall occurs. The oscillation frequency in 


Figure 4 is constant for the two cases, but Mach 
number and mean angle of attack are changed. The 
lift break occurs before the angle of attack reach- 
es its maximum value . In terms of the non-dimen- 
sional time parameter x*, the x* value of 6 at 
which lift stall occurs is reached before the maxi- 
mum angle of attack is reached. The Time Delay and 
Boeing Methods show comparable correlation for 
lift. For the Mach number 0.4 case (Figure 4b) the 
return to potential flow occurs earlier for de- 
creasing angle of attack than the return given by 
the Time Delay Method. Pitching moment correlation 
is better for the Time Delay Method. The triple 
loop characteristic is well duplicated. Similar 
correlation obtained with the Time Delay Method for 
a wide range of conditions demonstrated its promise 
as a practical method for analyzing unsteady aero- 
dynamics . 

Dynamic Stall Tests 

In order to obtain data useful in evaluating 
the two unsteady aerodynamic methods dynamic stall 
wind tunnel tests were run using a two-dimensional 
airfoil model. The model was oscillated at a 
frequency simulating the cyclic pitch variation on 
a helicopter rotor blade. Torsional frequencies 
representative of helicopter blade frequencies were 
obtained by varying a torsional stiffness element 
between the drive system and the airfoil section. 

The airfoil models were made to be as stiff as 
possible along their span and light in weight to 
approximate scaled helicopter blade mass and iner- 
tia properties. Hence the non-dimensional coeffi- 
cients in the equation of motion of the model air- 
foil were close to those of the helicopter blade 
torsional equation of motion based on the aero- 
dynamics of the three-quarter radius on the re- 
tracting blade. Two different airfoils were fab^ 
ricated, an NACA 0012 and an SC 1095. 

The model airfoils and drive system were 
designed to permit investigation of the effects on 
torsional response of torsional natural frequency , 
chordwise pitch axis location and torsional inertia 
over a range of IP frequencies for an NACA 0012 
airfoil and a cambered SC 1095 airfoil. The 
oscillating mechanism provided an 8-degree ampli- 
tude of motion of the model with an adjustable 
mean angle of attack. The model has a span of 
1.75 feet and a chord of 0.5 feet. The wind tunnel 
velocity was 275 fps for all tests. Time histories 
of the model angular motion were recorded at 
nominal driving frequencies of 8.0, 10.0, and 12.5 
cps . Tests were run for the full range of angle of 
attack for all the combinations of pitch axis, 
torsional inertia, torsional natural frequency, and 
airfoil type. A typical set of time histories 
for a basic reference condition (NACA 0012 airfoil, 
25 percent pivot axis, nominal blade inertia, and 
5P natural frequency ratio) is shown for four mean 
angles of attack in Figure 5- These time histories 
represent the time average of ten cycles . 

The elastic torsional deflection of the 
airfoil (difference between total angular motion 
and input angular motion) was obtained for each 
test condition by subtracting the input angular 
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Figure 5. Averaged Time Histories of Angle of 
Attack for the Model Airfoil. 

position time history from the averaged airfoil 
angular position time history: 


0(t) = ot(x) 


(a + TJsin27Ti) 

m 


where 0(x) is the difference between the measured 
non-dimensional angular time history response a(x) 
and the input driving motion. The non-dimensional 
time t is given by t/T, where the period, T, was 
established from the ten-cycle time-averaging 
process for that run. Some statistical variation 
in measured response was noted when stall flutter 
occurred, but in general the ten cycle time 
averaged response was representative of the 
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measured data. Two measures of stall response 
amplitude were extracted from each of the 0(x) time 
histories. These were A0-]_ which is one-half of 
the initial response to stall and 9 ijprp which is 
one-half of the overall peak-to-peak elastic de- 
flection. 

It was found that the initial stall response 
parameter A0, gave the most consistent indication 
of susceptibility to stall flutter. The possible 
reduction in flutter amplitude introduced by the 
time averaging procedure when there was cycle-to- 
cycle variation in phase made it somewhat difficult 
to assess the amplitude of flutter response. For- 
tunately, the initial stall deflection showed 
virtually no cycle-to-cycle variation. Figure 6 
compares measured initial deflection angles for an 
excitation frequency ft of 10 cps for the two air- 
foils at all combinations of airfoil natural 
frequency ratio (tj q = w Q /fti torsional inertia, and 
pitch axis. Certain general trends of deflection 
angle can be identified in the test results. 

1. Elastic deflection increases with mean 
incidence angle. 

2. For the same torsional inertia, response 
is generally greater for the lower 
frequency airfoil section. 

3. The amplitude of response is inversely 
related to torsional inertia. 

U. Forward movement of the pitch axis leads 
to a decrease in deflection. 

5. SC 1095 airfoil dynamic stall response 
begins to build up at a higher mean 
incidence angle than the 0012, but the 
two airfoils have comparable responses 
once stall is penetrated. 

Correlation Study of Two-Dimensional Results 

The two-dimens i on al flexured airfoil test 
data were compared with predictions based on 
various unsteady aerodynamic methods . The single 
torsional degree of freedom differential equation 
of motion for the flexibly mounted airfoil section 
oscillating in the wind tunnel test section is 
given by 

I-a.+ ca + Kta-Oju) = M(t) +K5' sin Sit 

where e = equivalent mechanical damping per unit 
span 

I = airfoil torsional inertia per unit span 
K = torsional spring constant 
Mfc= applied aerodynamic moment 
t = time 

« — airfoil angle of attack 
tx - amplitude of angular oscillation 
ow - mean angle of the oscillation 
ft - angular frequency of the applied torque 

This equation was solved numerically using 
the unsteady aerodynamic models to calculate the 
applied aerodynamic moment M(t). For the a. A, B 
Method unsteady data tables obtained from earlier 
oscillating airfoil tests. Reference U t were scaled 
for both airfoils. The measured steady state lift 
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and pitching moment data served as inputs in the 
Time Delay calculations. Additionally, the air- 
foil mean incidence angle used in the Time Delay 
solution was two degrees less than that set in 
the wind tunnel. The open jet flow deflection 
experienced at high unsteady lift coefficients was 
sufficient to decrease actual peak angles of 
attack to a value somewhat lower than the geo- 
metrically impressed pitch angle. The two-degree 
correction to gave consistently better corre- 
lation of the initial stall time. 

Correlation between measured wind tunnel 
model response and response calculated with the 
unsteady models was examined for thirty-six test 
conditions. The set of cases studied was suffi- 
cient to evaluate the independent effects on air- 
foil stall response of mean incidence angle. 


TEST 

,^.B 

TIME DELAY 




torsional natural frequency, chordwise pitch axis 
torsional inertia, and airfoil type. Relative to 
a baseline case taken to be the HACA 0012 airfoil 


at djj = lU°, cue /a = 5, TL. = 0.25, and I = I_, 
Figure 7 shows measured and predicted effects of 


mean angle, torsional natural frequency, pitch axis, 
and airfoil type on time histories of elastic de- 
flection. Coup ari son of the measured and predicted 
effects of airfoil mean angle of attack indicates 
that deeper penetration into stall results in 
sharper initial stall, deflection and larger resid- 
ual stall flutter response. The two analyses pre- 
dict these effects qualitatively, but each - 
especially the Time Delay model - overpredicts the 
amplitude of response. The main effects of an in- 
crease in torsional natural frequency are a shift 
in response frequency and a decrease in the am- 
plitude of elastic deflection. Figure 7 shows 
good correlation of response amplitude , although 
both analyses predict an initial stall response 
earlier than that measured. Moving the airfoil 
pitch axis forward causes delay in initial stall 
time and reduction in amplitude of response. The 
analytical results do predict the reduction in 
response amplitude, but the Time Delay model still 
results in overpredicted response. Finally, a 
comparison between the HACA 0012 and the SC 1095 
airfoils shows a delay in the initial stall time 
for the SC 1095 airfoil, which had a static stall 
angle measured in this wind tunnel to be about 
three degrees higher than that of the HACA 0012. 
However, the SC 1095 stall flutter amplitude was 
coup arable to that experienced by the HACA 0012 
at this condition. 


The time history correlation was good in 
that the initial response and the frequency of the 
subsequent oscillations were predicted. The trends 
observed in test were well matched by the analysis , 
although. the Time Delay model genera l ly overpre- 
dicted stall flutter response. The basic effects 
of structural changes on blade response time 
histories were well predicted hy either analysis. 


Although torsional elastic deflection is 
important in determining rotor stability and per- 
formance , the torsional moments resulting from 
stall flutter are the designer’s primary concern. 

To measure the trends of torsional moment with 
parameter changes, the twisting moment experienced 
by the flexible connector in the model airfoil 
drive system was calculated for each test condition. 
The torsion moment M@ was calculated using the 
equivalent spring stiffness of the connector: 

M e = V - (l airfoil “e 2) 9 


The torsional moments corresponding to the initial 
stall deflection angle were used to shew the 
effects of blade parameters on structural moments. 
Figure 8 presents typical results for three com- 
binations of airfoil type and mean angle of attack. 
It was generally found that decreasing torsional 
natural frequency reduced stall flutter moments. 
Although the stiffer system experienced lower 
response amplitudes, the corresponding structural 
moments were increased: 
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Figure 9 . Effect of Airfoil Type on Structural 
Moment s 


Forward placement of the airfoil pitch axis gener- 
ally decreased vibratory torsional moments. The 
two analyses predicted this trend with comparable 
accuracy. That forward movement of the airfoil 
pitch axis relative to the aerodynamic center 
reduces stall flutter moments can be understood 
based on lift and pitching moment hysteresis 
loops. For an airfoil with pitch axis forward 
of the center of pressure, positive lift forces 
cause negative moments about the pitch axis. 

For positive lift, the lift hysteresis loop is 
usually traversed in the clockwise direction, 
which contributes a negative pitching moment 
loop in the counterclockwise (stabilizing) 
direction. A decrease in torsional moment am- 
plitude with decreasing torsional inertia was 
generally found throughout the testing. This 
trend, evident in two of the conditions shown in 
Figure 8, is predicted somewhat more correctly 
by the Time Delay Analysis. Finally the two 
airfoils are compared in Figure 9« Fo** two 
different combinations of inertia and pitch axis, 
high stall flutter moments are delayed in mean 
angle with the SC 1095 airfoil. 


Flight Test Correlation 

Flight test data were correlated with the 
Normal Modes Blade Aeroelastic Analysis for both 
the CH-53A and CH-54B aircraft. Both models of 
unsteady aerodynamics were used. Information on 
the blade analysis used can be found in Refer- 
ence 10. 

Correlation of CH-53 control system loads, 
blade stresses and required power was studied at 
a nominal aircraft gross weight of 42, 000 lb 
(Cm/a = 0.083), a tip speed of 710 ft/sec, and a 
3000 ft density altitude for airspeeds ranging from 
100 knots to 170 knots. Inclusion of variable in- 
flow was found to be essential in calculating the 
proper levels of blade bending moments . It also 
provided some improvement in the correlation of 
blade torsional moments. 

The a. A, B and Time Delay aerodynamic models 
are compared at 137 knots in Figure 10. Figures 
10a and 10b shows that the computed blade stresses 
are comparable for the two methods. However, the 
push rod loads calculated with the Time Delay 
Model are much less than values calculated with the 
a, A, B Method and measured values. The Time 
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Figure 10. Correlation of CH-53A Blade Stresses and Pushrod Loads. 


Delay results generally do not give sufficiently 
large oscillations in stall. 

That better correlation of stall flutter 
moments was possible with the a. A, B Method is 
evident in Figure 11a which shows the buildup 
of vibratory pushrod load amplitude with airspeed. 
The a. A, B model predicts a buildup rate almost 
identical with the mean of the test data. A 
discrepancy of no more than 10 knots in the knee 
of the control load curve is evident at this thrust 
coefficient. Figure lib shows the correlation of 
pushrod load amplitude achieved with the a. A, B 
Method at three thrust coefficients. 

Calculated CH-5UB control loads were also 
generally less than measured values. Figure 12 
indicates that a definite stall boundary is pre- 


dicted by the analysis. Relative to the CH-53A 
calculations, a decreased control load stall 
speed and an increased rate of buildup with air- 
speed are clearly predicted. Again, higher loads 
are computed based on the a. A, B Model. The 
comparison of measured and predicted push rod 
load time histories indicates that the a. A, B 
results reflect a buildup in the higher fre- 
quency loads much more accurately than do the 
Time Delay calculations . 

It is not entirely clear why, relative to 
the a. A, B method, the Time Delay model under- 
predicts helicopter control loads while over- 
predicting the stall flutter oscillations of the 
two-dimensional wind tunnel model. Examination 
of several blade section pitching moment /angle 
of attack hysteresis loops indicates not so much 
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figure 11* Correlation of CH-53A Vibratory 
Pushrod Loads. 



Figure 12. Correlation of CH-5^B Vibratory 
Pushrod Lo ids. 


that more negative pitch damping is present in the 
ct, A, B results. Rather pitching moments along the 
blade are more in phase with each other, leading to 
larger modal excitation. In the a. A, B formu- 
lation, pitching moment coefficients are tabulated 
as functions of a, a and a values all along the 
blade. This formulation leads to similarly phased 
pitching moments. In the Time Delay Model, moments 
are calculated based on the angles of attack ex- 
ceeding the steady state stall angle for a certain 
interval of time and are not solely dependent on 
the instantaneous angle of attack characteristics. 
For small differences in calculated angles of 
attack, computed pitching moments for adjacent 
blade sections can be different in phase. Because 
the two-dimensional wind tunnel airfoil was modeled 
as a single panel for the calculation of aerodynamic 
forces, the effects of simultaneous spanwise stall- 
ing were not a factor in the correlation with that 
data. 

Comparison of Methods 

Because the a. A, B Method has demonstrated 
better correlation with flight test data, it con- 
tinues to be the method in use for blade design 
analysis. However, development of both methods 
continues. The a. A, B Method provides a relative- 
ly direct and simple procedure for calculating un- 
steady aerodynamic loads. Correlation has been 
good with test data but its disadvantage centers 
largely on the apparent need for extensive tests 
to provide the body of tabulated data required for 
each airfoil. Some success has been obtained by 
scaling the NACA 0012 unsteady aerodynamic tables 
based on steady state differences between airfoils. 
Work is also being done on developing analytical 
expressions to replace the tabulated data. These 
may lead to the ability to synthesize the data 
required for a given airfoil, which would make the 
method more desirable for general applications . 

The Time Delay Method has the great advantage 
of requiring only steady state airfoil data for its 
application. The correlation with forced oscil- 
lations of two-dimensional airfoils demonstrated 
its applicability over a wide range of conditions. 
Correlation with the tests described in this paper 
showed no clear advantage of the Time Delay Method 
over the a. A, B Method, and correlation with 
flight test data was definitely poorer with the 
Time Delay Method. Further work must be done to 
investigate the reasons for the poor flight test 
correlation. The problem may result from the 
assumption in the analysis that, on a blade, each 
radial section acts independently of its neighbor- 
ing sections. This causes a more random stalling 
along the span with time, which smoothes out the 
changes in blade loading. The propagation of stall 
along the span for the three-dimensional case of a 
helicopter blade must be added to the Time Delay 
Method. The ct, A, B Method does provide spanwise 
correlation in loading by use of torsional mode 
acceleration to calculate the B parameter. This 
acceleration is in phase for each point along the 
blade span. Incorporation of a suitable radial 
propagation model in the Time Delay Method may make 
this a more versatile, more easily applicable, and 
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more accurate model of unsteady aerodynamics. 

Until this can he shown the a, A, B Method 

continues in use in blade design. 

Conclusions 

1. The a. A, B and Time Belay unsteady aerodynamic 
models predict with good accuracy the lift and 
pitching moment hysteresis loops and the aero- 
dynamic pitch damping of rigidly driven os- 
cillating airfoils. 

2. Two-dimensional stall flutter tests indicate 
that reducing blade torsional stiffhess, re- 
ducing blade torsional inertia and moving 
blade pitch axis forward decrease stall flutter 
induced moments* Inception of stall flutter 
was delayed with the SC 1095 airfoil relative 
to the NACA 0012 airfoil; however, once 
initiated, stall flutter loads for the two 
airfoils were generally comparable. 

3. Stall flutter response of the two-dimensional 
model airfoils and the effects of airfoil 
structural design parameters on blade torsion- 
al moments can be calculated using both un- 
steady models. The Time Delay method gives a 
high prediction of response aaplitude. 

h. Good correlation of CH-53A and CH-5^B blade 
stresses and control loads was obtained with 
a rotor aeroelastic analysis employing vari- 
able rotor inflow and unsteady aerodynamics . 
Best correlation was achieved using the a, 

A, B unsteady model. The Time Delay method 
generally underpredicted full scale rotor 
stall flutter response . 

5. The a. A, B model is in use for blade design 
analysis. Refinements to the Time Delay Method 
may make it a more versatile and more easily 
applied unsteady aerodynamic model. 
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